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Abstract. Globular clusters produce orders of magnitude more millisecond pulsars per unit mass than the Galactic disk. 
Since the first cluster pulsar was uncovered twenty years ago, at least 138 have been identified - most of which are binary 
millisecond pulsars. Because of their origins involving stellar encounters, many of these systems are exotic objects that would 
never be observed in the Galactic disk. Examples include pulsar — main sequence binaries, extremely rapid rotators (including 
the current record holder), and millisecond pulsars in highly eccentric orbits. These systems are allowing new probes of the 
interstellar medium, the equation of state of material at supra-nuclear density, the mass distribution of neutron stars, and the 
dynamics of globular clusters. 
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INTRODUCTION 



BASIC PROPERTIES OF GC PULSARS 



The first globular cluster (GC) pulsar was identified 20 
years ago in the cluster M28 after intense efforts by an 
international team [1]. Since then at least 138 GC pul- 
sars', the vast majority of which are millisecond pulsars 
(MSPs), have been found. Finding these GC pulsars has 
required high-performance computing, sophisticated al- 
gorithms, state-of-the-art instrumentation, and deep ob- 
servations with some of the largest radio telescopes in the 
world, primarily Parkes, Arecibo, and the Green Bank 
Telescope (GBT). The payoff has been an extraordinar- 
ily wide variety of science. 

Low-Mass X-ray Binaries (LMXBs) have been known 
to be orders-of-magnitude more numerous per unit mass 
in GCs as compared to the Galactic disk since the mid- 
1970s [2, 3]. This overabundance is due to the produc- 
tion of compact binary systems containing primordially- 
produced neutron stars via stellar interactions within the 
high-density cluster cores. Since LMXBs are the progen- 
itors of MSPs, this dynamics-driven production mecha- 
nism also applies to them, and it has made GCs (par- 
ticularly the massive, dense, and nearby ones) lucrative 
targets for deep pulsar searches. 

Camilo & Rasio [4] produced an excellent review of 
the first 100 GC pulsars in 2005. This current review 
provides a significant update to Camilo & Rasio as it 
concentrates on the advances made (primarily with the 
GBT) within the past several years, including almost 40 
additional pulsars and over 50 new timing solutions. 



There are currently 138 known pulsars in 25 different 
GCs ^. Over 100 GC pulsars have been found in the past 
10 years, with almost 60 of these coming in the last 4 
years from searches using the GBT (see Figure 1). The 
three clusters Terzan 5, 47 Tucanae, and M28 account for 
approximately half of these pulsars, with 33, 23, and 1 1 
pulsars in each cluster respectively. Of the known pul- 
sars, 80 are members of binary systems^, 50 are isolated, 
and 8 are as yet undetermined. 

Almost 90% of GC pulsars are true MSPs with spin 
periods Pspi„ < 20 ms (see Figure 2). Their spin proper- 
ties seem consistent with the standard "recycling" sce- 
nario [e.g. 8], with surface magnetic field strengths of 
B ^ 10^^^ G and characteristic ages of Tc ^ lO^^'^yrs. 
Most of the rest of the pulsars are partially recycled. 
However, there are several seemingly very out-of-place 
"normal" radio pulsars with P(p,„ > 0.2 s and 10^ yrs 
as well [9]. 

There are at least four distinct groups of binary GC 
pulsars (see Figure 3). The first two are similar to the bi- 
nary MSPs found in the Galactic disk. The "Black Wid- 
ows" have very low mass companions (Mc ^ 0.04 M©) 
and orbital periods of several hours, while the "nor- 
mal" low-mass binary MSPs (LMBPs) likely have He- 
lium white dwarf (WD) companions of mass Mc ^ 0.1 — 
0.2 Mq and orbital periods of several to tens of days. The 
Black Widow systems are relatively much more common 



' For an up-to-date catalog of known GC pulsars, see Paulo Freire's 
website at http : / / www . naic . edu/~pf reire/GCpsr . html 



^ The Galaxy has roughly 150 known GCs [5]. 

' One of the "binaries", B 1620— 26 aka M4A, is in fact a confirmed 
triple system including a white dwarf and a planet [6, 7]. 



138 pulsars in 25 clusters 



CO 

o 



O 

cu 

_Q 



o 



o 



o 



O 



O 
1985 



Ln 
oo 



1^ o 



Telescope PSRs 



Jodrell Bank 
Parkes 
Arecibo 
GBT 
GMRT 



1990 



(N 

CD 
O 
O 



1^ 

CD 
O 



ro 

in 

CD 

o 
o 



o 

ro 



5 

40 
27 

65 




0°^ 



9 



1995 
Year 



2000 



in 1^ 

CD CD 

o o 

CD CD 



O 
CD 
1^ 
CD 
O 
CD 



CN 

in 

CD 
CJ 



CD 
O 
CD 



in 

CD 
O 
CD 



O 2 

fli 



I 



1^ 



2005 



cn 

CM 



CN 

in 

CD 
O 
CD 



CD 



O 

CD 
O 



CD 
CD 
O 
O 




I 




o 

ro 



O 



Globular Cluster 

FIGURE 1. The number of pulsars per globular cluster and ( inset) a timeline of globular cluster pulsar discoveries. 



in GCs {^25% of the binaries) than in the Galactic disk, 
though (~4% of the binaries). 

The other two groups of binary GC pulsars are possi- 
bly unique to clusters and their formation therefore likely 
depends on the high stellar densities and interactions 
found in the cores of GCs. Approximately 10% of GC 
binaries appear to have "main sequence"-like compan- 
ions which show irregular eclipses, erratic timing, and 
often have hard X-ray and/or optical counterparts. The 
prototype system is J1740-5340 in NGC 6397 [10]. Fi- 
nally, ~20% of the known GC binaries have highly ec- 
centric orbits (with e > 0.1). The standard recycling sce- 
nario produces circular orbits due to tidal interactions 
during mass transfer, and so the large eccentricities are 
probably either induced during multiple stellar interac- 
tions with passing stars [e.g. 11] or produced directly 
during an exchange encounter with another star or bi- 
nary. Over the past several years the numbers of pulsars 
in each of these two groups have grown dramatically. 
This is likely due to the fact that recent surveys have suc- 
cessfully probed many of the most massive and dense 
clusters in the Galaxy where these systems are preferen- 
tially produced (i.e. M28, M62, 47 Tucanae, Terzan 5, 



NGC 6440, and NGC 6441). 



SEARCHING FOR CLUSTER PULSARS 

The signal-to-noise ratio from a radio obser- 
vation towards a (perhaps unknown) pulsar is 
o<: 5vAe v^^Tfofei- Sy is the pulsar's flux density 
at the observing frequency v. Ag is the effective collect- 
ing area of the telescope. r„.s is the system temperature, 
which at the ~GHz radio frequencies of interest is 
roughly the sum of the receiver temperature (typically 
15—25 K), the 2.7 K cosmic microwave background, and 
the Galactic synchrotron background (7bo/o=v~^-^ and 
typically a few to tens of K depending on sky position 
and observing frequency v). By is the radio bandwidth 
used for the observation and t^bs is the observation 
duration. In addition, signal-to-noise ratios improve 
when a pulsar has short duration pulsations compared 
to its pulse period (i.e. a small pulse duty cycle). If 
the pulses are smeared or broadened in time, perhaps 
by uncorrected orbital motion or interstellar medium 
effects, signal-to-noise ratios during searches can be 



reduced to effectively zero . 
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FIGURE 2. Spin frequency histogram of the 138 currently 
known GC pulsars. Eighty of the pulsars are confirmed mem- 
bers of binaries, 50 are isolated, and 8 are as yet undetermined. 
It is interesting to note that the binary MSPs seem to spin more 
rapidly on average than the isolated MSPs. Perhaps this is an 
indication that they are in general younger (i.e. more recently 
recycled) than the isolated MSPs. Such an explanation makes 
sense if all isolated MSPs originally come from binaries and 
therefore must destroy their companions over time. 
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FIGURE 3. Orbital period Pg^i, vs. minimum companion 
mass {Mc^min'i assuming a pulsar mass of I.4M0) for the 70 
binary GC pulsars with well-determined orbits. All 6 pulsars 
with "main sequence"-like companions and 14 of the known 
pulsars in highly eccentric (e > 0.1) orbits are labeled, as well 
as many more "normal" binary systems. Pulsars in "NGC" 
clusters are labeled without "NGC" to save space. The large 
grouping of pulsars with Pgyi, < 1 day and Mc.min < 0.04 Mg 
are the so-called "Black Widow" systems. The lack of pulsars 
in the lower right portion of the diagram is not due to selection 
effects, as those pulsars (if they existed) would be relatively 
easy to identify during searches or through timing observations. 
M15C is the only confirmed GC double neutron star system 
[12] while TerSN is likely the only known GC pulsar with a 
Carbon-Oxygen WD companion [13]. 



The Problems with Cluster Distances 

JVISPs are intrinsically very faint radio sources. Be- 
cause of this, all of the wide-area Galactic pulsar sur- 
veys conducted to date have been severely sensitivity 
limited for MSPs (many have been instrumentation hm- 
ited as well). Of the ~60 known Galactic TvISPs, ~80% 
are within 2 kpc of the Sun. In contrast, the nearest GCs 
(MA and NGC 6397) are just over 2 kpc from the Sun, 
and most, including the best targets for pulsar searches, 
are at distances of 5 — 15 kpc [5]. 

In addition to the inverse square law problem, the large 
distances to GCs often (especially for the bulge clusters 
near the Galactic center) imply large column densities of 
the interstellar medium (ISM). The ionized ISM causes 
frequency-dependent dispersion of radio waves v~^), 
scatter-broadening of the radio pulses (o= v '* '*), and for 
certain clusters dramatic fluctuations of observed pulse 
intensity due to diffractive scintillation [e.g. 15]. These 
ISM effects, as well as a substantially reduced Galac- 
tic synchrotron background, have pushed typical observ- 
ing frequencies from ^400 MHz in the early 1990s up 
to 1.3— 2GHz in the past decade. At these frequencies, 
especially with much wider observing bandwidths avail- 
able (hundreds of MHz), significant sensitivity gains 
have been realized despite the usually steep radio spec- 
tra of the pulsars themselves (flux densities Sy oc v" with 
-3^ a <.-l and <a>~-1.8; [16]). 



Searching for Binaries 

Since most GC pulsars are in binaries^, orbital mo- 
tion causes Doppler variations of the observed pulsa- 
tion frequencies during an observation. If uncorrected, 
these variations can make even very bright MSPs unde- 
tectable. However, correcting for unknown orbital mo- 
tion identically would be extremely computationally ex- 
pensive and so current searches only account for linear 
changes in apparent spin frequencies (i.e. constant /). 
These "acceleration" searches [e.g. 17] are valid when 
the orbital period is much longer than the observation 
duration {Porb^^Qtobs)- 

Acceleration searches add an extra dimension to the 
traditionally two-dimensional phase space of dispersion 



* See Lorimer & Kramer [14] for a much more thorough discussion of 
pulsar search sensitivities. 

' The current binary fraction of ~60% is a lower limit since finding 
isolated pulsars is much easier than finding binaries. 



measure DM^ and spin frequency / over which one 
must search for pulsars. Since the number of trials in 
the acceleration or / dimension is proportional to f^^j, 
the long observations used to improve the sensitivity of 
GC searches greatly increase their computational costs. 
Typically, searches for the first pulsar in a cluster are 
made using a large range of likely DMs but only a limited 
range of accelerations or /. Once the first pulsar is found 
and the rough DM toward the cluster is known, a much 
smaller range of DMs is searched, but with a much larger 
range of possible accelerations for additional pulsars. 

As an example, to properly search a single 7 hr GBT 
observation of Terzan 5 (where the DM is known to 
~5%) with a fuU range of acceleration searches, requires 
approximately one CPU-year of processing on state-of- 
the-art CPUs. However, it is important to realize that 
without acceleration searches (or other advanced binary 
search techniques such as Dynamic Power Spectra [18]), 
the majority of the binary GC MSPs that have been 
uncovered over the past decade would simply not have 
been found. 



A Renaissance in 2000 

The above paragraphs summarize why GC pulsar 
searches require long integrations (sensitivity) at GHz 
frequencies (minimize ISM effects), using the largest 
telescopes (collecting area), the best receivers (wide 
bandwidths and low system noise), and large amounts 
of high-performance computing (acceleration searches). 
The dramatic increases in the numbers of known GC 
pulsars begirming in 2000 (see Figure 1) resulted from 
significant improvements in each of these areas. First, 
new low-noise and wide-bandwidth (Bv~300MHz) ob- 
serving systems centered near 1 .4 GHz became available 
at Parkes and Arecibo. Second, the rise of affordable 
cluster-computing allowed acceleration searches to be 
conducted at investigator institutions rather than at spe- 
cial supercomputing sites. The first major success from 
these improvements (and a significant driver for further 
efforts) was the discovery of 9 new binary MSPs in 
47 Tuc L15J. 

The third and perhaps most important improvement 
was the completion of the GBT in 2001. With its state- 
of-the-art receivers, approximately three times greater 
than Parkes, and the ability to observe over 80% of the 
celestial sphere, it is perfectly suited to make deep ob- 
servations of GCs. By the end of 2003, a fantastic wide- 
bandwidth (Bv'~600MHz) system centered near 2 GHz 



* DM is the integrated electron column density along the Une-of-sight 
to a pulsar. 



became available which provided 5—20 times more sen- 
sitivity for MSPs in certain GCs in the Galactic bulge 
than the 1 .4 GHz system used at Parkes. The discovery of 
30 new MSPs in Terzan 5 [13], a cluster previously ex- 
tensively searched at Parkes, including the fastest known 
MSP (J1748-2446ad aka TerSad; [19]), were some of 
the first results. Pulsar surveys of many additional GCs 
using the same system are ongoing. 



Future Cluster Pulsar Surveys 

Recent work on the luminosities L of GC pulsars [20] 
has confirmed earlier results [e.g. 21] suggesting that 
the luminosity distribution roughly follows a d log = 
—d log L relation. In addition, this work implies that we 
currently observe only the most luminous pulsars in each 
cluster. Together, these facts indicate that our current GC 
pulsar surveys are completely sensitivity limited such 
that even marginal improvements in search sensitivities 
will result in new pulsars^. The history of GC pulsar 
searches has directly demonstrated this fact many times. 

However, there seems to be httle hkehhood of making 
very large improvements in GC pulsar search sensitiv- 
ities over the next several years. Most of the variables 
in the signal-to-noise equation are already nearly opti- 
mal (e.g. Tsys, By, tobs^ v). Dramatlc improvements 
in sensitivities and therefore pulsar numbers will almost 
certainly require a new generation of larger telescopes 
(i.e. larger Ae) such as FAST^ or the SKA'. 

WHICH CLUSTERS HAVE PULSARS? 

Figure 1 shows the number of pulsars in each of the GCs 
with known pulsars. Currently there are 10 clusters with 
5 or more pulsars and 3 clusters with 10 or more pul- 
sars: M28 with 11, 47 Tucanae with 23, and Terzan 5 
with 33. Camilo & Rasio [4] pointed out that there are 
very few clear correlations between cluster parameters 
and the numbers of known pulsars. In fact, the only sim- 
ple properties that seem to be related to the number of 
known pulsars are the total mass of the cluster (which 
likely influenced how many neutron stars were originally 
retained) and the distance D to the cluster (since all GC 
pulsar searches are currently sensitivity limited). How- 
ever, even these indicators have exceptions. For example. 



' 47 Tuc is a possible exception to this rule as deep radio imaging [22] 
and X-ray observations [23] indicate that scintillation may have aheady 
allowed the identification of nearly all of the observable MSPs in the 
cluster. 

FAST: http ; / /www . bao .ac.cn/LT/ 
' SKA: http://www.skatelescope.org 



CO Centauri has been searched extensively but unsuccess- 
fully with the Parkes telescope, yet it is one of the nearest 
and most massive GCs in the Galactic system. 

A more sophisticated indicator of which clusters may 
contain more LMXBs and therefore MSPs is the pre- 
dicted stellar interaction rate in the cores of the clus- 
ters, where LMXBs and MSPs are likely formed. Pooley 
et al. [24] showed a strong correlation between the num- 
ber of X-ray sources in a cluster and its Tc, which they 
expressed as Tc °= Pg ^r^, where po is the central density 
and rc is the core radius. We can attempt to adjust the 
indicator to account for our senstivity issues by ranking 
clusters by TcD^^. Using this metric, we find that many 
of the clusters with numerous pulsars are near the top 
of the list, including 47 Tuc, Terzan 5, M62, NGC 6440, 
NGC 6441, NGC 6544, M28, and M15. Also near the top 
are several others which hkely contain numerous pulsars 
but whose positions near or behind the Galactic center re- 
gion (and therefore large amounts of ISM) make searches 
very difficult (e.g. NGC 6388 and Liller 1 [25]). 

It is important to reaUze, though, that because of lim- 
ited amounts of telescope time, pulsar searchers have 
specifically targeted those clusters near the top of the 
YcD^ list first. Therefore, clusters further down that list 
may simply not have known pulsars because they haven't 
been searched to the same sensitivity levels as the clus- 
ters near the top of the list. Until we have a large number 
of GCs, independent of their position on the TcD^^ list, 
searched to sensitivities comparable to the recent GET 
2 GHz surveys, it will be difficult to determine just how 
good of a predictor TcD '^ really is. 

The recent 1 .4 GHz survey of all 22 GCs within 50 kpc 
and visible with Arecibo [20] is a good example of the 
type of surveys we need. The Arecibo survey found 1 1 
new MSPs, the majority of which are in clusters with 
fairly average values of r^D^^. No pulsars were found 
(or have ever been found) in clusters with very low cen- 
tral luminosity densities, po < 10^L©pc^^. A similar 
unbiased survey of ^^60 GCs at 1 .4 GHz using the Parkes 
telescope has uncovered 12 new pulsars [10, 26] in 6 
GCs. Unfortunately, the limited sensitivity of that survey 
does not rule out even relatively bright MSPs in many 
of the clusters, thereby making it difficult to draw con- 
clusions about cluster properties and their pulsar popula- 
tions. 



What pulsars are in those dusters? 

We can compare the pulsar populations in the best 
studied clusters, such as Terzan 5 and 47 Tucanae, to at- 
tempt to determine if the properties of the clusters affect 
their pulsars. Two of the simplest things to compare are 
the spin-period distributions and the binary populations. 



both of which do show significant differences. 

The spin-period distribution of the Terzan 5 pulsars is 
significantly flatter than that of the 47 Tuc pulsars with 
more faster and more slower pulsars. In fact, Ter 5 con- 
tains 5 of the 10 fastest known spinning pulsars [19], 
while 47 Tuc contains only one within the top ten. Like- 
wise, 47 Tuc has no MSPs rotating slower than ~8 ms, 
whereas Ter 5 has six. A Kolmogorov-Smirnov test sug- 
gests a <10% chance that the two period distributions 
were drawn from the same parent distribution. As for the 
binary populations, Ter 5 has only two known "Black- 
Widow" systems compared to five in 47 Tuc. On the 
other hand, Ter 5 has six highly-eccentric binaries com- 
pared to none in 47 Tuc. Are these differences related to 
the current interaction rates in the cluster cores (Ter 5's 
is 2—3 times that of 47 Tuc's) or perhaps the "epochs" 
when the MSP production rates were the highest? 

The pulsars in M28, NGC 6440, and NGC 6441 are 
more similar to those in Ter 5 than to those in 47 Tuc, 
both in terms of their spin periods (i.e. flatter distribu- 
tions) and their binary parameters (with a broader mix 
of different systems). It is interesting to note, though, 
that the 10 known pulsars in the very similar clusters 
NGC 6440 and NGC 6441 rotate on average significantly 
slower than those in either Ter 5 or 47 Tuc. Only one pul- 
sar rotates faster than 5 ms, and half rotate slower than 
13 ms [27]. 



XmiNG OF CLUSTER PULSARS 

While it is obviously essential to find the pulsars in GCs 
to do any science with them, without detailed follow-up 
observations, and in particular pulsar timing solutions, 
the amount of science one can do is severely Umited. For 
most GC pulsars the extraordinary precision provided by 
MSP timing provides ~0.1" and often significantly bet- 
ter astrometric positions (crucial for multi-wavelength 
follow-up [e.g. 23]), extremely precise Keplerian or- 
bital parameters for binaries, and measurements of the 
pulsar's apparent spin period derivative. For some pul- 
sars, particularly those in eccentric orbits, certain post- 
Keplerian orbital parameters can be measured [e.g. 14] 
which allow the determination of (or at least constraints 
on) the masses of the pulsar and/or the companion star. 
Establishing timing solutions for as many GC pulsars as 
possible allows us to use them both individually and in 
ensembles to probe a wide variety of both basic physics 
and astrophysical phenomena. Over the past two years, 
the number of GC pulsars with timing solutions (cur- 
rently 107) has almost doubled, resulting in many inter- 
esting new results. 




FIGURE 4. Positions for tlie 32 Terzan 5 pulsars witli timing 
solutions overlaid on 1.4 GHz radio contours from the VLA 
[25], and an NTT I-band image (courtesy S. Ortolani; [31]). 
The average positional error in RA is <0.02" while for DEC it 
is <0.2" (the much larger DEC errors are due to the fact that 
Terzan 5 is very near the ecliptic). The good correspondence 
between the timing positions and the VLA radio contours im- 
plies that potentially all of the radio flux is indeed produced 
by unresolved radio pulsars (many tens of which are as yet 
unknown). 

Ensembles of pulsars 

The largest ensemble of pulsars in a single GC with 
timing solutions are the 32 in Terzan 5 (only TerSU, a 
weak eccentric binary, remains without a solution; see 
Figure 4). However, there are five other GCs with at least 
five pulsars with timing solutions as well (47 Tuc [28], 
M28 [29], 1V115 [21], NGC 6440 [27], and NGC 6752 
[30]). These ensembles of pulsars in individual GCs can 
produce unique science. 

Probes of Ionized ISM and Intra-Cluster Medium. 
The precise DMs (errors <0.1pccm^^) and timing po- 
sitions for 32 pulsars in Ter 5 have recently allowed a 
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FIGURE 5. Radial distribution comparison of "isolated" pul- 
sars (including binaries with very low-mass companions, < 
0.05 Mp ) and heavier binary pulsar systems. The numbers on 
the thin black lines are various q = Mp/M^ values where Mp 
is the mass of the pulsar system and is the mass of the 
dominant stellar component in the GCs {M^, ~ 0.9 M r,). Sur- 
prisingly, the "isolated" systems appear to be more centrally 
condensed (and therefore possibly more massive) than the bi- 
nary systems. 



unique probe of the ionized Galactic ISM between us and 
the cluster on parsec scales [32]. A calculation of the DM 
structure function indicates that the fluctuations in the 
ISM on 0.2— 2 pc scales roughly follow those predicted 
for Kolmogorov turbulence. Earlier work on 47 Tuc us- 
ing the pulsar positions, accelerations (see Figure 6), and 
DMs provided the first definitive measurement of ionized 
gas within a GC [33]. 

Statistical Neutron Star Mass Measurement. Using 
the projected offsets of pulsars from their cluster cen- 
ters and a model for how relaxed components of different 
masses should be distributed within a cluster, it is possi- 
ble to statistically measure the masses of the pulsar sys- 
tems Mp [34]. Figure 5 shows the 107 GC pulsars with 
timing positions split into two different groups: 1) iso- 
lated pulsars or binary pulsars with very low-mass com- 
panions (Mcjnin <^ 0.05 M©) and 2) binary pulsars with 
more massive companions. Surprisingly, the "isolated" 
systems seem to be more centrally condensed than the 
(supposedly more massive) binary systems. Fits of the 
observed distributions give q = Mp/Mj, ~ 1.5 for the bi- 
naries and ^ ^ 1.7 for the "isolated" systems. Assuming 
that the dominant stellar components in the cluster cores 
have mass = 0.9 M0 implies binary system masses of 
Mp^hin ^ 1.35 M0 and larger masses Mpj^o ^ 1.53M0 
for the "isolated" pulsars. The reason for this difference 
in mass segregation is currently unknown. 



Cluster Proper Motions. The very precise positions 
available from MSP timing allow individual pulsar 
proper motions given regular observations over 5—10 
years. The measurement of several pulsar proper motions 
from a single cluster allows a measurement of the proper 
motion of the GC itself. Currently this has been accom- 
plished for three clusters (47 Tuc [28], M15 [12], and 
NGC 6752 [30]), and several more will likely be mea- 
sured within the next couple of years. Such measure- 
ments are very important for determining the Galactic 
orbits of GCs and predicting the effects of tidal stripping 
and/or destruction. Measuring cluster proper motions is 
very difficult in the optical (using HST, for instance), es- 
pecially for the Galactic bulge clusters which are distant 
and plagued by extinction. 

Cluster Dynamics. The projected positions of the 
pulsars with respect to the cluster centers as well as 
measurements of their period derivatives (which are usu- 
ally dominated by acceleration within the gravitational 
potential of the GC; see Figure 6) provide a sensitive 
probe into the dynamics of the cluster and even con- 
strain the mass-to-light ratio near the cluster center [35]. 
These measurements can provide evidence for the pres- 
ence (or absence) of black holes in the cores of the clus- 
ters [36] . Cluster dynamics also influences pulsars by 
ejecting some of them to the outskirts of the clusters 
or even entirely [37]. A recently uncovered example of 
such a system is M28F, a bright (for a GC MSP) iso- 
lated pulsar located almost 3' from the center of M28. 
That offset is larger than for any other GC pulsar except 
for NGC 6752A [30], where exotic ejection mechanisms 
have been invoked to explain its position [e.g. 38]. 



Individual Exotic Pulsars 



'Young' Pulsars with t 10'' yrs 



, , n n n , i H , , n 

10 ' 2x10 ' 3x10 ' 4x10 ^ 

Pulsar 'Acceleration' = cP/P (m/s^) 

FIGURE 6. Histogram of observed globular cluster pulsar 
"accelerations" {cPohsl Pobs)- For GC pulsars, the observed ac- 
celeration is the sum of the pulsar's intrinsic acceleration (i.e. 
cPint/Pint), an apparent acceleration due to the pulsar's proper 
motion, and accelerations from both the Galaxy's and the glob- 
ular cluster's gravitational potentials [35, 39]. The proper mo- 
tion and Galactic potential terms are typically small compared 
to the others. Since GC pulsars are typically fully recycled with 
small intrinsic accelerations (~ 10~' m s~-^), the GC gravita- 
tional accelerations (~ |10"^|ms"^) usually dominate. This 
is apparent in the figure by a nearly symmetric clustering of 
pulsars around zero acceleration. The pulsars on the observer's 
side of a cluster receive positive accelerations while those on 
the far side receive negative accelerations and appear to spin 
more rapidly with time. The five pulsars with observed accel- 
erations > 10^^ ms^^ are anomalously young (characteristic 
ages 1 X 10^ < Tc < 3 X 10^ years) and their intrinsic accelera- 
tions are much larger than the maximum possible gravitational 
acceleration from their clusters. Apparently clusters continue 
to produce such systems [e.g. 37]. The three pulsars with the 
most negative accelerations provide unique probes of the cen- 
tral dynamics (and lower limits on the mass-to-light ratio) of 
their parent clusters [e.g. 21]. 



While ensemble studies of GC pulsars are very inter- 
esting, many of the pulsars are truly exotic or unique ob- 
jects, worth studying individually. 

Young and Slow Pulsars. A handful of slow, "nor- 
mal" (Tc ~ lO' yrs) pulsars have been known in GCs 
for some time: B1718-19 aka NGC 6342A [40], 
B1820-30B aka NGC 6624B [41], and B 1745-20 aka 
NGC 6440A [9]. Recently, at least one more slow pulsar 
has been uncovered, NGC 6624C with P^pm = 0.405 s 
[18], as well as M28D, an 80 ms binary pulsar that is 
definitely "young" (see Figure 6). The slow pulsars, 
which have likely not been through any recycling, must 
have formed relatively recently even though all of the 
massive stars would have gone supernova lO'^'yrs ago. 
One possibiUty is that the pulsars formed via electron- 
capture supemovae, perhaps via accretion-induced 
collapse of a massive WD or merger-induced collapse of 



coalescing double WDs [see N. Ivanova's paper in these 
proceedings; 37]. 

The fastest MSPs. Ter5ad is the fastest MSP known, 
with a spin period of Pspm = 1.396ms [19]. Its discovery 
finally broke the 23-yr-old "speed" record estabUshed by 
the very first MSP discovered [42] and renewed hope for 
finding a sub-MSP (a pulsar having a spin period under 
1 ms). A sub-MSP would provide by far the most direct 
and interesting constraints on the properties of matter 
at nuclear densities and would be of major significance 
to physics in general [43]. Besides TerSad, several other 
very rapid rotators have been uncovered in Ter5 (Ter50 
at 1.676ms and Ter5P at 1.728ms; [13]), and it seems 
likely that the first sub-MSP (if they exist) might be 
found in a GC. 



"Main-Sequence"— MSP Systems. Several pul- 
sars have been recently discovered (including 47TucW 
[44, 45], M62B [46], Ter5P [13], TerSad [19], andM28H 
[29]) which appear to have bloated "main-sequence"- 
like companion stars much like the prototype system 
J1740-5340 in NGC 6397 [10, 47]. These pulsars are 
eclipsed for large fractions of their orbits and show 
irregular eclipses on some occasions. Timing positions 
usually associate them with hard X-ray point sources 
where the high-energy emission is likely generated via 
colliding MSP and companion winds. Several of the 
pulsars have been identified in the optical where they 
exhibit variability at the orbital period. In addition, 
at least some of them exhibit highly erratic orbital 
variability (resulting in several large amplitude orbital 
period derivatives) likely due to tidal interactions with 
the bloated companion stars. These systems could be 
the result of an exchange encounter between a main 
sequence star and a "normal" binary MSP system. 
Alternatively, perhaps the companions are the stars that 
have recently recycled the pulsars and we are observing 
newly born MSPs [10, 47]. Multi-wavelength studies of 
these systems are difficult [e.g. 45], but allow a wide 
variety of additional constraints to be placed on MSP 
emission mechanisms, their winds, and the evolutionary 
histories of their systems (including tidal circularization 
theory). 

Highly Eccentric Binaries. At least 15 GC pulsar 
systems are members of eccentric binaries with e > 0.1, 
and most of those contain MSPs. In contrast, only a sin- 
gle eccentric binary MSP is known in the Galactic disk'". 
Ten of these systems have been discovered since 2004: 
six in Terzan 5 [13], two in M28 [29], and one each in 
NGC 1851 [48, 49] and NGC 6440 [27]. Given the an- 
gular reference that ellipses provide, pulsar timing can 
easily measure the precession of the angle of periastron 
or d). For compact companions, d) will be dominated 
by general relativistic effects, and its measurement pro- 
vides the total mass of the binary system [14]. Timing 
observations of four of these systems (Ter5I & J [13], 
NGC 6440B [27], and M5B [Freire et al. in prep.]) indi- 
cate "massive" neutron stars (>1.7M0) which constrain 
the equation-of-state of matter at nuclear densities [43]. 
Such constraints are impossible to achieve in nuclear 
physics laboratories here on Earth. In addition, similar 
measurements for M28C, a 4.15 ms pulsar in an 8-day 
orbit with e = 0.85, indicate that the pulsar is less mas- 
sive than 1.37M0. This is a fairly low mass for a neu- 
tron star which must have accreted a substantial amount 
of material during recycling and will likely constrain re- 



See D. Champion's contribution to these proceedings. 



cychng models. For more information on these systems, 
see P. Freire's contribution to these proceedings. 

Other Exotica. There is already one confirmed GC 
triple system (PSR B1620-26 in M4) which contains 
an MSP, a white dwarf, and a planetary-mass compo- 
nent [6, 7]. Intriguingly, ongoing observations show very 
strange and seemingly systematic timing residuals from 
the "isolated" MSP NGC 6440C [29]. One explana- 
tion for these residuals is the presence of one or more 
terrestrial-mass planets. Given the strange variety of sys- 
tems that have already been found in GCs, it is quite pos- 
sible that one of the many currently uncharacterized sys- 
tems could be another unique object. 

PROSPECTS FOR THE FUTURE 

Given the wide variety of science (most of which was 
unanticipated) that has already resulted from GC pulsars, 
and the fact that we are currently only seeing a small frac- 
tion (perhaps 10-20%) of the total GC pulsar population, 
the future of GC pulsar astronomy seems very bright. 
Improvements in search sensitivities with current instru- 
ments will likely uncover tens of additional systems 
and next-generation telescopes like FAST or the SKA 
promise to find hundreds. Among these pulsars may be 
even more spectacular "exotica" such as MSP— MSP or 
MSP— black hole binaries, and from these new pulsars 
will come many surprising results. 
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